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Abstract 



Ultraluminous X-ray sources fed by Wolf-Rayet star winds and X-ray bursters in ul- 
tracompact binaries with He or C white dwarfs have accretion disks, whose properties 
I may significantly differ from those of pure H a-accretion disks. We have therefore in- 

ly-^ ■ eluded the dependence on c harge number Z and mean molecular weights ^^/i into the 

I IShakura and SunvaevI l)l973j) scaling relations for the key parameter s of the disk. Further- 

Q^ ■ more, we also consider the case of the pseudo- Newtonian potential of lPaczvhskv and Wiital 

0^ I (|l980|). These scaling relations might become useful, e.g., when estimating the illumination 

i eflSciency of the external parts of the disk. We also address the changes in the structure 

I of the boundary (spreading) layer on the surface of neutron stars, occurring in the case of 

. H depleted accretion disks. 

^ ■ 
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! 1 Introduction 
> . 

^ , Recent observations of galaxies by the CHANDRA X-ray observatory confirm the existence of 

^ I compact, ultraluminous X-ray sources (ULXs), emitting at luminosities considerably higher 

than the Eddington luminosities of neutrons stars or stellar mass black holes with 1 — 15 Mq 



( Boberts et al. . l2fln4al lbl: ILui and Mirabel 120051 ). The most drama tic results w ere obtained 



by CHANDRA observations of star forming galaxies and regions (,Grimm et al.l . 12003.) , where 
over the large luminosity range L = 10^^ — lO^'^ erg/sec the corresponding high mass X- 
ray binary (HMXB) luminosity function is well fitted by a single power law with a strong 
cut-off at L ~ lO^*^'^ erg/sec, which is more than two orders of magnitude higher than the 
Eddington luminosity of a solar mass neutron star. Due to the fact that the HMXB luminosity 
function does not exhibit a change of slope or other peculiarities in the region just above the 
Eddington limit of stellar mass black holes, it is plausible to assume that ULXs represent 
the high-mass/high-accretion-rate tail of the 'ordinary' HMXB population with masses in the 
range 10 - 15 Mq (King, 2002; Grimm et al. 2002; Gilfanov, 2004). 

These new observational results increased strongly the theoretical interest in the ques- 
tion, whether an accreting object can produce super-Eddington luminosities, or whether the 
Eddington luminosity provides a real upper limit for the luminosity of accreting objects. In 
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principle, one can think of several different mechanisms leading to super-Eddington luminosi- 
ties and, at present, there does not yet exist a commonly accepted model for the huge energy 
output of ULXs. For example, young rotation powered pulsars or jet sources directed towards 
us could subs tantially increase the obs erved X-ray flux. Other possibilities might be slim ac- 
cretion disks (jAbramowicz et al.','l988'), or an accretion disk that experiences radiation driven 
inhomogeneities (Begelmann . 2002)- Alternatively, many observers and theorists believe that 
the high lumin ositv of ULXs is due to accretion onto an int ermediate mass black hole with 
102 - 10^ Mq dMiller and Colbertl 120041 : iKubota et aP . l2002l 'l. 

On the other hand, there also exist several 'minor' effects (e.g., related to a disk's in- 
clination and chemical composition), each of which being, in principle, able to increase th e 
observable luminosity of an accreting object by a factor of two or three (jGrimm et al, ]. l2002l V 
Due to the aforementioned homogeneity of the HMXB luminosity function, it is worthwhile to 
carefully reconsider the contributions of these 'minor' effects, even though most of them are 
known to theorists for more than 30 years. In the present paper we are going to focus on the 
case, where the chemical composition of the accreting matter is strongly deviating from the 
standard cosmic abundances. Obviously, in this case the mass per electron in the accreting 
gas can be higher than for a pure H plasma or standard cosmic abundances. It is well-known 
that for completely ionized He, C, O, N, or Mg plasmas, having twice less electrons per baryon 
than H, the Eddington luminosity is by a factor of two higher than for a pure H plasma. 

Prom the observational point of view there are evidences that, under cer tain conditions . 



accreting gas in binary systems can consist only of elements heavier than H ()Hammer et al 



hm. For- example, several observations of X-ray sources ^n extremely dense Unary systems 
reveale d a neutron star accreting matter supplied by a He, C, O, Ne , or e ven Mg white 
dwarf ( Schulz et al. , 200ll : Juett and Chakrabartv , 200,4 Nelemans et al. , 20041 ) . The unusual 
chemical abundance influences the power release due to nuclear explosions on the surface of 
the neutron star, which are observable as X-ray bursts. The duration of bursts and their 
repetition rate are also affected strongly. 

The problem of non-standard abundances might also be of relevance, if the accretion onto 
the black hole of the HMXB takes place in a star forming region or galaxy. Then many of 
the HMXBs are fed by stellar wind accretion. A particular interesting scenario in this regard 
corresponds to the case, where the donor star is of Wolf-Rayet-type. During their evolution 
process these stars may have lost a significant amount of their H rich envelope, leading to 
stellar winds domina ted by He or, in some case, even bv C and N 

2000l : Ivan der Huchti . I2OOII : ICrowther ei~al1 . 12OO2I : lAbbott et all . .2004 ). For instance, one of 
the brightest, most luminous X-ray sources in the Milky Way, Cygnus X-3, is fed by the dense, 
high-velocity wind from its He Wolf-Rayet companion (van Kerkwijk et al., 1992; Lommen et 
al., 2005). 

In the case of accretion onto a black hole, the plasma density in the region of the 
main energy release of a standard H dominated accretion disk is not high enough to pro- 
duce the photons necessarv for creating a black body radiation spectrum inside the disk 
( Shakura and Sunvaev . 1973h . Moreover, in the limit M/Me 1 and a — > 1 the optical 



depth becomes rather small (here, as usual, Me denotes the accretion rate corresponding to 
the Eddington luminosity and a is the viscosity parameter). It is, therefore, the purpose of 
the present paper to study how the disk's density, temperature, and optical depth depend on 
the chemical abundance in the extreme case of H poor accretion. By simply repeating the ap- 
proach of Shakura and Sunyaev (1973), we will show that under such peculiar conditions the 
density of electrons and nuclei increases, thereby also increasing the effective optical depth. 
This might stabilize the disk against the transition to the two-temperature hot flow regime 
( Shapiro et"Zl.ll97(j'l. 

Recentlv. iHammerlHIl B have performed numerical simulations to determine the 
structure in the outer regions of C/O/Ne-dominated accretion disks in ultracompact binary 
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systems (for details of their NLTE-model see Nagel et al., 2004). Their primary objective was 
to explain the observed peculiarities in the optical spectra of these systems, by taking into 
account metal-line blanketing and irradiation by the central object. Ho wever, in spite of con - 
siderable theoretical interest in the properties of H poor accretion disk (jMenou et al. 



we were unable to find in the literature simple analytical results, describing the power law 
scaling relations for key parameters of the standard accretion disk with respect to its chemi- 
cal composition in the extreme case of H or even He depleted accretion flows. Therefore, the 
formulae given below, which have originally been calculated as part of a more complex study, 
will possibly be useful when interpreting observational data and deriving simple analytical 
estimates. Finally, it is also worthwhile to mention that a high abundance of He or heavier 
elements does not only decrease the disk (scale) height by a factor of about two, but also 
changes the width and the surface brightness of the spreading layer on the surface of neutron 
stars. Therefore, in the periods between X-ray bursts, the accretion of heavier elements also 
influences the heating of the outer parts of the disk by X-ray radiation emitted from the 
central disk regions (this al so applies to the case of acc retion onto black holes) and from the 
surface of the neutron star ( Lvutvi and Sunvaev , 197(i ) . 



2 Structure of helium or metal dominated accretion disks 



The hydrodynamic equations considered by IShakura and Sunvaev do not explicitly 



depend on the chemical composition of the disk. More precisely, in their model the chemical 
composition enters only via the thermodynamic equations for the pressure and the energy 
density of the gas: 

_ gkT _ 3 gkT 

finip 2 finip 

and via the Rosseland mean opacities for Thomson and free-free processes: 

CTT 0.398 2 _i 
KT = = cm^ g , (lb) 

= „,2^fi^V?V"'-=g- (10) 



/Xe/i/ \ nip J \K 

Here Z denotes the charge number of the ions, p the mass density of matter, k the Boltzmann 
constant, cjt the Thomson scattering cross section, T the temperature, nip the proton mass, 
/ie/7 the mean molecular weight of the electrons and ions, respectively, and = + fij^ . 
The results for the disk's scale height, surface density, etc., g iven below in Eqs. 



were obtained by using Eqs. and following the approach of lShakura and Sunvaev! (|l973l ^. 
In particular, Eqs. @ - @ refer to a Newtonian potential $ = —GM/{E?' + z^)^/^, where 
R is the cylindrical radius. It is however easy to see that the scaling relations with respect 
to the chemical parameters (see also Table remain also valid, when including relativ istic 
corrections in the form of the pseudo-Newtonian potential of lPaczvnskv and Wiital ()l98n'V as 



evident from the results in Appendix A. It is convenient to express the results in terms of 
non-dimensional quantities m = M/Mq, rri = M /M^, r = R/tq and s = 1 — r~^/^, where 
ro = 3rG = QCM/c? = ni 8.86 km. The total luminosity of an infinite a-disk, ranging from 
ro to oo reads Lqo = C-^c^ with an efficiency = 1/12 in the case of a Newtonian potential 
and Q ~ 0.06 for a Schwarzschild black hole. The critical Eddington luminosity, providing an 
upper limit for spherically symmetric accretion, can then be written as 

2TTrG ni„c^ . _i X , s 

Le = Aie^E,H = fJ-e — = fJ-em (1.26 X lO^"* erg sec ^), (2) 
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where, here and below, quantities labelled by H refer to the case of a fully ionized, pure H 
plasma. The corresponding critical accretion rate, obtained from the condition Lqo = L^, 
reads 

Me = MeME,H = /^e m (2.22 x 10"^ Mq year"^). (3) 

For the inner disk zone (A), characterized by dominating radiation pressure Pg <C P-y and 
^ Kg, one finds 

H = mrhs (2.2 km) (4) 

T. = C fie a^^ m'^ (9.9 x 10^ g cm'^) 

Qc = a'^ m'^ m'^ r^/^ s"^ (2.2 x 10"^ g cm"^) 

riec = m'^ r^/^ s"^ (1.3 x 10^° cm"^) 

VR = a r~^/^ s (-2.5 x 10^ km sec~^) 

e^c = 7^-3/2 ^^^^ ^ -j^qIS ^^.^ cm~^) 

Pc = A^e r~^/2 (1.0 x 10^^ erg cm"^) 

T, = ^-1/4 ^-3/8 (2.5 X 10^ K) 

^-1/16 ^-2 ^93/32 ^-2 4 ^ iq-2)^ 

where we note that rh = mn/ fi^- Here the index c is used to label midplane values, and H is 
the disk scale height, S the surface mass density, vr the (vertically averaged) radial velocity, 
e^c the energy density of radiation, Uec = Qc /il^enT-p) the number d ensity of electrons, and 
r = (^T^ff)^^^ S/2 the effective optical depth (jShakura and Sunvaev . 1973h . 
Similarly, for the intermediate zone (B) with Pg ^ P^ and ^ Kg: 

H = r'/V"'/vy'%~'/'° (5) 



m 



9/10 



^1/5 ^21/20 ^1/5 ^ 10-1 j^^^) 



m^/^ ^3/5 ^-3/5 ^3/5 7 ^ g ^^-2^ 
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1/5 ™2/5 -2/5 -3/5 



s"^/^ (-5.3 km sec"^) 



m ' m ' r 
^-4/5 ^8/5 ^-18/5 g8/5 ^ ^ jgl^ erg cm" 



= C-"/5 ^2/5 ^9/10 ^-9/10 

^-9/10 ^4/5 ^-51/20 ^4/5 (2.3 x lO^^ erg cm"^) 

^-1/5 ^2/5 ^-9/10 ^2/5 (12 X 10^ K) 
Te = Z r ^/1%21/10 ^-1/5 ^^-1/2 ^-4/5 
1/5 ^1/10 J/20 1/10 



m ' m 
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Table 1: Relative quantitative changes compared to a pure, fully ionized H accretion disk 
with [ii = [If. = Z = \ and same parameter values (a,M, M). Inner zone A: Pg <C -By and 
KT ^ /tfi- Intermediate zone B: Pg ^ P-y and kt ^ Kfr- Outer zone C: Pg ^ P-y and 
kt ^ /tfT- All values refer to pure disk plasmas, containing o nly a single ion species. Th e 
same results are obtained for the pseudo-Newtonian potential of Paczviiskv and Wiital (jl98nl ). 
see Appendix A. 
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and for the outer zone (C) with Pg ^ Py and ^ Kff: 
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The disk zones are labelled by A, B, and C as in Shakura and Sunyaev (1973). Numerical 
values based on the above results for different types of pure, fully ionized plasmas are given in 
Tabled For example, this table indicates that, compared to a pure H disk with same param- 
eter values (a, M, M), the disk becomes typically thinner by a fa ctor of two, if dominated b y 
heavier elements. As one may check by virtue of Saha' s equation ( Landau and Lifshitd . 20031 ) . 
considering fully ionized He, C, Mg, etc. is a good approximation in the case of ultracompact 
binaries possessing hot accretion disks with small radial diameters. 

In addition to the quantities from Eqs. (jlj) - B , we listed in Table |T] the relative changes 
of the characteristic time scale for Comptonization ( Pozdnvakov et al. . Tl983h 
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Another important time scale is giv en by the equ ipartition time for the energy exchange 
between fast electrons and slow ions (jSpitzei . E3) 
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(8) 



where n/c = Qc/{l^imp) is the number density of ions in the disk midplane. The quantities 
teq and tco can be compared with, e.g., the time scale of radial motions = R/vr, or, 
alternatively, with the revolution time scale tn = 2ttR/Q. It is worthwhile to note that, in 
contrast to ieq/co/R) the revolution time scale ta is independent from the chemical composition. 
A two-temperature flow regime can devel op in the hot, in ner disk region, if the latter becomes 
optically thin and if teq ^ tQ holds ((Shapiro et al. , 19761 ) . 



6 



3 Boundary (spreading) layer of neutron stars 



The main difference between accretion onto black holes and neutron stars is due to the fact 
that on the neutron star's surface the kinetic energy of the accreting matter is transformed 
into radiation, whereas in the former case this energy is absorbed by the black hole. There 
exist several different models for the boundary layers of neutron stars (Sunyaev and Shakura, 
1986; Popham and Sunyaev, 2001). In the spreading layer model of Inoeamov and SunvaevI 
(iiii), ;wo bright belts, located equidistant from the equator on the surface of the (slowly 
rotating) star, emit radiation corresponding to the local Eddington flux 



qEie) 



GM nipc 



(9) 



where i?* denotes the radius of the neutron star, 6 the latitude angle (0 = corresponds 
to the 'equatorial' disk mid-plane), vk = (GM/R^,)^/"^ is Keplerian velocity and f,^ is the 
rotation velocity on the surface of the star (to be obtained by solving the boundary layer 
equations of Inogamov and Sunyaev, 1999). The local Eddington flux qEid) is determined by 
the difference of the gravitational force and the centrifugal force due to the rotation of the 
accreting matter on the surface of star. Hence, because of Le = fieLE,H, formula © implies 
that the flux that may be emitted from the surface of a neutron star increases with /ig (small 
deviations from this simple proportionality may occur due to a weak dependence of v^f, on the 
chemical parameters and /i). Therefore, in the case of He or metal rich accretion flows, 
this flux is increased by a factor of two (at same values of M), compared to a pure H disk. 
Moreover, similar to the scale height H of the disk, the meridional hei ght -ffsT, of the boundar 



(spre ading) layer, which can be estimated from the energy balance (llnogamov and Sunvae 



ES) 



iry 



LsL = Le 



/ie-^E,H 



ii* 



(10) 



decreases at same values of M, if the disk is dominated by elements heavier than H. 
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A Paszynsky-Wiita potential 

As discussed by IPaczvnskv and Wiital (|l98(]l 'l , the effects of general relativity can be modeled 
by replacing the Newtonian potential $ = —GM/r with the modified potential 

GM , , 

a>i?,z = , 11 

r — tq 

corresponding to an efficiency factor = 1/16, in contrast to the Newtonian value C, = 1/12. In 
order to simplify subsequent formulae, it is convenient to define three dimensionless correction 
factors by 

R r , , 

Xo = ^ = -j::, 12a 

R — rc r — 1/3 

_ 3 R + ^/rn^ _ 3 r + ^ 

^^-2 2{R-rG) 2 2(r-l/3)' ^ ^ 

« - {-^-^^- 3-) - (- ^ -) ■ (-) 

As can be seen from these factors, relativistic corrections are important only in the innermost 
parts of the disk. Therefore, we will confine ourselves here to giving the results for zones A 
and B. 



Zone A: In the case of dominating radiation pressure, Pg ^ -By, and dominating Thomson 
scattering, kt Kfj, one finds 



H{R) 


= Xo^ X2 H 


t(R) 


2 -2 
= Xo Xl X2 ^ 


Qc{R) 


= Xo Xl X2 ^ 6c 


hec{R) 


= Xo Xl X2 ^ '^ec 


vr{R) 


= Xo ^ Xi^ xl VR 


e-yc{R) 


= xlxiX'.'e,^ 


Pc{R) 


= xlxix-,'Pc 


fc{R) 


1/2 1/4 -1/4 
= Xo Xl X2 


fc(R) 


25/8 17/16 — 
= Xo Xl X2 



where on the rhs. the results from Eqs. @ are to be inserted (here and below non-hatted 
quantities refer to the expressions derived for the case of the Newtonian potential). 
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Zone B: In the case of dominating gas pressure, Pg ^ P^, and dominating Thomson 
scattering, kt ^ Kfr, one finds 

(14) 





-1 1/10 1/10 

— Xo Xl X2 


u 
n 




4/5 -1/5 ^ 
— Xl X2 ^ 






7/10 -3/10 

— Xo Xl X2 


8c 




7/10 -3/10 
— Xo Xl X2 


flee 




-4/5 1/5 - 

= Xl X2 VR 






4/5 4/5 

= Xl X2 ^ic 




Pc{R) 


9/10 -1/10 
= Xo Xl X2 


Pc 


Tc{R) 


— Xl X2 




fc{R) 


1/2 4/5 -7/10 
= Xo Xl X2 ^c, 



where on the rhs. one has to insert the results from Eqs. ©• 

A dimensionless transition radius r^B, separating zones A and B, can be determined from 
the condition f^{rAB) = f^{rAB) ; i.e., TAB is given by the solution of the equation 

r21/16>,-1^5/4 ^1/8 ^^-9/8 ^ 3 ^_i ^1/2 ^-1/8 ^1/8 ^1/8 (^5^ 

The function on the Ihs. attains its minimum value 14.8 at r = 2.54. Consequently, zone A 
exists only if 

0.02 < /xV2 ^-1/8 a^/s ^1/8 (Ig) 
where we have already inserted the efficiency factor C, = 1/16. 
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